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ABSTRACT 


Offshore jets or filaments are a well observed phenomenon in the California 
Current System. Their surface and subsurface structures as well as their spatial and 
temporal variability are well documented. The methods by which these filaments are 
formed is not as well understood as their structural form. The most common 
generation theories are: (1) variations in wind stress coupled with topographic 
irregularities, (2) dynamic instability and (3) geostrophic turbulence. In this study an 
attempt is made to identify the factor(s) responsible for filament formation south of 
Cape Mendocino by numerically simulating the California Current System using a two- 
layer, non-linear, primitive equation model. It is shown that baroclinic instability is 
the primary method by which filaments are formed in this region and that other factors 
such as barotropic instability, friction and non-linearities can alter the characteristics of 


the fastest growing baroclinically unstable wave. 
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1. INTRODUCTION 


A. THE CALIFORNIA CURRENT SYSTEM 

The Cahfornia Current System, CCS, has been extensively studied by Elickey 
(1979). The CCS consists of four separate currents that are spatially and seasonally 
MistmbUted Gi etiem, este Gast ol the Wimited States. They are: (1) the California 
Current, a southward flowing surface jet; (2) the California Undercurrent, a northward 
flowing sub-surface jet; (3) the Davidson Current, a poleward flowing surface jet north 
of Point Conception; and (4) the Southern California Countercurrent, a northward 
flowing surface current south of Point Conception in the California Bight (Hickey, 
19779). Figure I.) shows the average seasonal configurations from [Ikeda and Emery 
(1984) for the CCS. Typical surface jet velocities are on the order of 20 to 40 cm/s. 
The California Undercurrent has a tvpical strength of about 5-10 cm/’s. For the 
profiles on the left, where the California Undercurrent is present, the surface jet depth 
is about 150 m while for the profiles without a subsurface jet, the surface jet extends 


down to about 300 m. 


EARLY SUMMER LATE SUMMER / WINTER + SPRING 
iN Oo 20cemea"! s FALL 





Figure 1.1 Seasonal Configurations of the California Current 
System (Ikeda and Emery, 1984). 


B. FILAMENTS IN THE CALIFORNIA CURRENT SYSTEM 
The presence of offshore jets or filaments in the CCS is well documented. 


Satellite radiometry and color scanning along with intensive ia situ data such as that 


obtained during the Coastal Ocean Dynamics Experiment (CODE) have provided large 
data bases for these features. 

From an observational standpoint, filaments are well documented. Definitive 
Studies by Flament er al. (1985), Kosro (1987) and Kosro and Huyer (1986), for 
example, have greatly expanded our knowledge of filament surface and subsurface 
Structures along with their spatial and temporal variability. 

Filaments consist of strong baroclinic jets flowing offshore at peak velocities of 
0.5 m’s or greater (Rienecker er al., 1985), with narrow return flows to the south 
(Flament ez al., 1985). The offshore extent of filaments is typically on the order of 200 
to 300 km with maxima of around 500 km (Brink and Hartwig, 1985). Filaments 
advect cold upwelled coastal waters away from the continental shelf. 

Filament widths are on the order of 50 km (Flament ef al., 1985). The 
wavelength associated with their occurrence at numerous points along the coast is 
roughlv 100 to 400 km for the region to the north of Cape Mendocino (Ikeda and 
Emery, 1984). Formation time scales are generally on the order of weeks with lifetimes 
on the order of a month. Figure 1.2 shows the surface velocity field derived from 
satellite feature tracking for a filament studied by Flament ez al. (1985). 

Filaments are seasonal in nature. Few have been observed during late winter 
through early spring. Their first appearance is off the Baja Peninsula in March. 
Between April and May they occur between Point Conception and Cape Mendocino 
and by late summer they reach their maximum northward extent off the Oregon coast. 
Generally, filaments have been more observed south of Cape Mendocino. (Brink and 
Hartwig, 1985) 


C. FILAMENT FORMATION 

The precise formation process for filaments is not as well understood as their 
characteristics once formed. Numerous theories abound with the most prevalent 
causative forces being: (1) coastally linked processes associated with variations in wind 
stress coupled with irregularities in the alongshore topography, (2) dynamical 
instabilities of the CCS, and (3) geostrophic mesoscale eddy field interactions with the 
CCS. It is possible and indeed probable that more than one theory may be responsible 
for filament formation. 

Ikeda and Emery (1984) demonstrated by a combination of observational data 
and numerical experiments that coastal irregularities off the coast of Northern 


California and Oregon are strongly correlated with filament formation. They also 
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Figure 1.2 Surface Velocity Field of a Filament (Flament er al., 1985). 


showed that baroclinic instability was probably the primary causative force north of 
Cape Mendocino. Batteen (1987) demonstrated that baroclinic instability could occur 
in the CCS due to imposing an unstable jet profile similar to the early summer profile 
shown in Figure 1.1. Chen (1986) showed that baroclinic instability coupled with 
topographic variations are important to the destabilization of the CCS in the vicinity 
of Cape Mendocino. Mysak (1977) also demonstrated the importance of baroclinic 
instability in a numerical study of the CCS off Vancouver Island. Mooers and 
Robinson (1984) documented an eddy pair off Northern California, while Owen (1980) 
and Huyer et al. (1984) observed strong eddy-like features in the northeastern Pacific 


close to the continental slope. 


D. PURPOSE OF THIS STUDY 

The objective of this study is to determine the causative factor(s) involved with 
filament formation along the coast of California south of Cape Mendocino. A two- 
layer, non-linear, primitive equation model with topography is used to simulate the 
dynamics of the CCS. The results strongly point to baroclinic instability as the 


primary generation mechanism for filaments in this region. 


1] 


The specifics of the model and the various numerical experiments are covered in 
Chapter II. The linear baroclinic instability problem is formulated in Chapter III and 
the results applied to the data from the numerical experiments to verify the claim of 
baroclinic instability to the CCS. In Chapter IV, the results are discussed and 
conclusions presented. 
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If. NUMERICAL SIMULATION OF THE CALIFORNIA CURRENT 


A. THE NUMERICAL MODEL 
1. Model Equations 

Simulation of the CCS is performed using a two-layer, primitive cquation, 
Semu-implicit numerical scheme initially developed by Flurlburt (1974). This scheme 
has been employed in numerous ocean mesoscale circulation studies (e.g., [lurlburt and 
i@tomposon, 1980) 1952; Smuth and O Brieng 1983; and Smith, 1986). Linear test cases 
have been run for comparison with linear analytic solutions to show that this approach 
1s valid (Smith and Reid, 1982). Motion in each laver is governed by a momentum 
equation (2.1) and a continuity equation (2.2). 


av. 
Wig (VeVitVieV)vy + KX £Vy = -hy VP; + A,VVy 2.1) 


vt 


3h. 
<i, gey 
Ce 


I 
O 
~_ 
bo 
tJ 
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Where 1 = layer index (1= 1 upper, 1= 2 lower) 


Variables and notation are defined in the Appendix. The fluid is assumed to be 
metrocstatic and Boussinesq. [he density in each layer is constant. Ihe effects of 
winds, tides, thermodynamics and thermohaline nuxing are not included. 

2. Nlodel Domain 

A rectangular region (1100 km X 800 km) was divided mnto 10 km by 10 km 
squares to form the grid for the numerical model finite differencing. By rotating the 
grid through a desired angle, the model can be used to simulate the flow within the 
CCS as indicated by Figure 2.1. 

Topography is included by approximating the major features apparent in the 
Mendocino Escarpment region. Figure 2.2 shows the digitized topography in the 
vicinity of the Mendocino Escarpment (Chen, 1986) and Figure 2.3 shows the 
topography used in the model. Principal features of the topography used in the model 


are: (1) a shelf region of 500 m depth extending out 50 km from the coastline which ts 
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Figure 2.1 Geographical Positioning of Model Domain. 


Solid box represents a basin rotation (y) of 90°, dashed box for rotation of 120° 
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a uniform vertical wall; (2) a shelf slope region extending from 50 km to ~ 100 km with 
an exponentially decreasing slope; (3) a modeled Mendocino Escarpment (without a 
ridge) which separates a 3000 m depth region to the north from 4500 m water depth to 
the south. [he topography rotates with the model domain. 
3. Boundary Conditions 

Using the orientation of Figure 2.1, no-slip boundary conditions are set on the 
eastern and western boundaries. The northern boundary condition consists of a 
prescribed inflow in the upper layer. For simulated lower layer flows to the south or 
north the lower layer boundary consists of either a prescribed inflow or a radiation 
condition (Camerlengo and O’Brien, 1980) respectively. The southern boundary 
condition is set up just opposite the northern. The upper layer always has a radiating 
condition and the lower is radiating or prescribed for southward and northward lower 
layer flows respectively. 

Obviously, these boundary conditions do not reflect totally the situation in the 
CCS. By putting a closed boundary to the west, there will inevitably be shear 
produced vorticity that is not realistic. However, it is believed that the east-west 
dimension of the domain is large enough so that the vorticity production at the western 
boundary will not affect the dynamics of the CCS which is located near the eastern 


boundary. 


B. NUMERICAL EXPERIMENTS 
1. The Reference State 

In order to filter through the many possible generating mechanisms for 
filaments, variation of selected parameters is necessary. The problem is twofold: (1) 
selective variation of parameters may cause model deviations from the CCS 
configuration such that extrapolation of results from the model to the CCS is 
impossible. The results would still be important from a dynamical standpoint but 
Without the connection with the ‘real world’ their usefulness is limited; (2) isolation of 
specific model responses by’ changes in model parameter(s) can be inaccurate due to 
complex cross correlations between varied and non-varied parameters. 

In order to assess the model output variations resulting from parameter 
changes, a reference state is established. Instead of a simple configuration such as flat 
bottom barotropic flow, the reference state is selected to be as complex as possible 
involving such things as baroclinicity, topography and the presence of a pre-existing 


eddy (all suggested to be dynamically important in Chapter 1). As mentioned before, 
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pography in Vicinity of Mendocino Escarpment (Chen, 1986). 
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Smoothed Topography used in Numerical Simulations 


applicability of the results to the CCS is vital and use of a reference state as near as 
possible to the CCS configuration keeps the study focus on CCS dynamics rather than 
general dynamics. 

The reference state consists of a two-layer geostrophically balanced, Gaussian 


jet of the form, 
hi(y)=A.{l-exp(-y?/2L7)] i= 1,2 (2.3) 


where L is the e-folding width scale of the jet and the y-axis oriented east-west as 
shown in Figure 2.1. The flow strength is determined by the amplitude coefficient A; 
which determines the magnitude of the surface and interface height anomalies through 
hy). Additionally, a two-layer geostrophically balanced, axisymmetric, Gaussian eddy 


is included, 
hi(x,y)=A-{l-exp(-R?/2L*)] i= 1,2 (2.4) 


where R*=(x*+y?) and represents the radial distribution of the eddy. 

Mean upper layer thickness is set at 150 m (H,) which corresponds to a first 
internal Rossby radius, Ry, of ~18 km. This is comparable temthe saa Rg for the 
CCS (Brink and Hartwig, 1985). H, is dependent on the topography depth, D(x,y), by, 


H,(x,y) = 4500 - H, - D(x, y) (22) 


The jet e-folding scale, L, is set at 25 km. The jet is fixed parallel to the coast over the 
slope region and the eddy positioned offshore from the jet. The spatial relationships 
between the jet, eddy and topography are indicated in Figure 2.4. In Figure 2.4 and in 
all subsequent figures of model output fields, tick marks along the x and y axis 
represent a distance of 275 km and 200 km respectively. The lines of constant surface 
height anomaly can be thought of as streamlines representing the barotropic part of 
the current. The interface anomaly is sloped opposite the surface with an amplitude 
such that its baroclinic component in the opposite direction sets up the selected flow 
magnitude in the lower layer. 

The initial jet and eddy strength is set at 0.4 m/s (southward for jet and 


anticyclonic for the eddy) in the upper layer and 0.01 m/s in the same direction in the 
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Figure 2.4 Reference State Flows and Topography. 


Reference state initial conditions of (a) surface height een and $b) height, of 
Pemom aoove 4500 m (m}. Contour intervals are (a) 2.5 cm. and (b) 250 mn. Lick 
marks represent 275 km in the x direction and 200 km in the y direction. 
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lower layer. The horizontal Laplacian frictional coefficient, Ay, is set at 500 m?/s 
which is roughly in the middle of the range of A, (100-1000 m*/s) typically used in 
numerical studies. A summary of model parameters for the reference state is shown in 
Table I. 


TABEE | 
REFERENCE STATE MODEL PARAMETERS 


PARAMETER SYMBOL VALUE 
E-W Basin Extent be 800 km 
N-S Basin Extent be 1100 km 

| Initial Upper Layer Thickness H, 150 m 
Initial Lower Layer Thickness H,(x,y) Variable 
Maximum Basin Depth (Hy + Ho) ax 4500 m 
Depth of Bottom D(x,y) Variable 
Coriolis Parameter (@ 40° N lat.) fy 0.94 x 10457} 
dfidy Bo 20 x 107! mls 
Gravitational Acceleration g 98ms7 
Reduced Gravity 9’ 2.0 x 10° ms? 
Horizontal Eddy Viscosity Coef. Jay 5.0 x 102 m* s7! 
Tinee:S-ep At 3600 s 
Grid Increment in x-direction Ax 10 km 
Grid Increment in y-direction Ay 10 km 
Jet Maximum, Upper Layer Vij 0.4ms! (southward) 
Jet Maximum, Lower Layer V3; 0.01 ms! (southward) 
Eddy Maximum, Upper Layer Vie 0.4ms"! (anticyclonic) 
Eddy Maximum, Lower Layer Vow 0.4ms! (anticyclonic) 
Jet Position Offshore Lj 50 km (from east wall) 
Eddy Position Offshore ve 200 km (from east wall) 
Basin Rotation Angle y 90° (from east) 


Figures 2.5 and 2.6 show the reference state model output fields for height 


anomalies and upper layer velocities respectively. For this and all subsequent model 


output field plots, surface height anomaly is contoured in cm, interface height anomaly 
is contoured in m and the velocities are contoured in cm/s. Based on the directional 
characteristics and the magnitudes of the flows, it is apparent that some phenomenon 
very similar to filament formation is occurring. Of particular note is the velocity 
structure of the offshore flow versus the inshore flow. The offshore flow has strong 
unidirectional flow while the inshore flow is slower and has more directional variability. 
The maximum offshore velocity by day 30 is about 50 cm/s and the inshore flow 
maximum about 40 cm/s. This type of structure and magnitude of flow correlates well 
With observed data in filament studies (Flament er al., 1985; Kosro, 1987). 

There are also two important signatures of baroclinic instability in the output 
fields: (1) wave growth in the initially nearly quiescent lower layer is indicative of 
vertical energy transfer from the upper to lower layer, and (2) there is a 180° phase 
shift between the waves in the upper layer and the waves in the lower layer. 

2. Variation of Parameters 

Structures shown in Figure 2.5 are consistent with observations of the form 
and magnitude of filaments. In order to isolate potential causes for filaments, the 
assumptions that the wave structures in Figure 2.5 represent filament formation and 
growth are made. Since the reference state is as complex as possible, the variations of 
the numerical simulations are largely one of simplification in specific areas of the model 
to see whether filaments are altered or made non-existent. 

All of the numerical simulations are for 30 days and certain parameters from 
Table 1 are never varied. Table 2 provides a list of those parameters which are varied 
in the different experiments and what their reference state value is. 

a. Experiment No. I (H, Increased to 250 m) 

Upper layer thickness was increased to 250 m to determine the model 
Sensitivity to the vertical structure. Figure 2.7 shows the model output fields for 
surface and interface height anomalies. From the fields, it 1s apparent that the wave 
Structures are significantly different than those in the reference state. They are weaker 
and seem to stop growing halfway through the sequence. This points to the 
importance of the vertical structure to the wave development process. 

b. Experiment No. 2 (Eddy Removed) 

The eddy was removed in this experiment to determine whether it was 
responsible for the wave structures by its interaction with the coastal jet. Figure 2.8 


shows the output fields for surface and interface height anomalies. The results are very 
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Figure 2.5 Model Output Fields of Height Anomalies 
for Reference State: 


Height anomaly fields for (a) surface at day 15 om) (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.6 Model Output Fields of Upper Layer Velocity 
for Weference State. 


Velocity fields (cm/s) for (a) upper layer at day 15 and (b) upper layer at day 30. 


zy 


TAD bee 
VARIATION OF MODEL PARAMETERS 


PARAMETER VALUES (Notation from Table 1) 


ble exG yy) V4; V9; Vie) or L; ee Y 
Exp. # 
rel 150 ] -4 = -.0] -4 -0O1 50 200 =90 
I 250 l -4 +0] -4 -0O1 50 200 890 
2 150 l -.4 -.O1 ne ne “Woo 200 =90 
3 150 l -4 -.01 -4 = -.01 180 =6300-—Ss( 90 
4 150 0 -4 = -.01 -4 -01 50 200 90 
5 150 ] -4 = -,05 -4 -05 50 200 90 
6 150 ] -.4 Ol -.4 Ole se 200 890 
7 150 ] -4 -.01 -4 -01 50 200 0 
8 150 ] -4 = -.0] -.4 -O1 50 200 ~=180 
o 150 0 -4 -01 ne ne 50 200 90 
10 150 ] -4 +01 ne ie wo0 200 =120 
11 150 | -.2 -.01 -2 -O1 50 200 90 
12 150 ] -.] -.O] -.1 -O1 50 200 90 
Notes: Depths in m 

Lengths in km 
| Velocities in m/s, + to north, - to south 


Eddy rotation: CW for + vel., CCW for - vel.; ne= no eddy 


| 

2 

3 

: Bottom depth: 1 topography 0 flat bottom 
3 Basin rotation angle (y) in degrees from east 


similar to the reference state. The wave structure 1s more uniform in this run, 
indicating that the eddy is not causative. However, it does affect the wave structure. 
c. Experiment No. 3 (Jet Moved Farther Offshore) 
The jet and the eddy were moved out farther away from the coast to see if 
Wave structures were being generated as a result of boundary interactions with the 
eastern, no-slip boundary. Figure 2.9 shows the output fields for this run. In general, 


the results are very comparable to the reference state. It looks as though the full 
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Figure 2.7 Experiment # | (H, Increased to 250 m) 
Model Output Height Anomalies. 


Fleight anomaly fields for (a) surface at dav 15 yom (b) surface at day 30 (cm), (c) 
interface at dav 15 (m), and (d) interface at day 30 (m). 
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Figure 2.8 Experiment # 2 (Eddy Removed) 
Model Output Height Anomalies. 
Height anomaly fields for (a) surface at da 


é 15 yam), (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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growth and extent of the waves to the west is now being hindered but their structure is 
consistent. Apparently, the boundaries can restrict but cannot generate the 
disturbances. 

d. Experiment No. 4 (Topography Removed) 

In this case, topography was removed and a flat bottom at 4500 m depth 
utilized. There is substantial circumstantial evidence linking filaments to coastal 
topographic features such as the Mendocino Escarpment. By removing the 
topography, it is hoped that it can be determined whether bottom variations are 
causing the filament formation. Figure 2.10 shows the output fields for this run. The 
fields are nearly identical with the reference run. Apparently, topography 1s not 
inducing the filament formation process for this weak lower layer flow. 

e. Experiment No. 5 (Increased Lower Layer Velocity) 

In this run the lower layer velocity was increased from | cm/s to 5 cm/s 
flowing southward, which is considered to be on the order of the maximum observed 
velocity in the deep water south of Cape Mendocino (Stabeno and Smith, 1987). This 
reduces the vertical shear but increases significantly the model coupling with the 
bottom. Although the previous run tended to discount the importance of the bottom, 
it is important to see what possible effects it might have in a worst case situation. 
Figure 2.11 shows the output fields over the domain. The results are nearly identical 
with the reference state with the amplitudes slightly smaller. This indicates that the 
model is relatively insensitive to southward lower layer velocity. 

f. Experiment No. 6 (Lower Layer Flow Direction Reversed) 

In this experiment, the lower layer flow was reversed so that it was flowing 
Opposite the upper layer flow. This effectively increases the vertical velocity shear and 
examines a possible configuration of the CCS where there is a deep, weak northward 
flow. Figure 2.12 shows the output fields for this run. The wave structure is similar 
to, yet noticeably stronger than, the reference. This would indicate the wave growth to 
be very sensitive to the vertical velocity shear. 

g. Experiment No. 7 (Zonal Flow Westward) 

In order to perform a sort of check on the model dynamics, the basin was 
not rotated into its CCS position. Since most instability studies are concerned with 
zonal rather than meridional flows, the analytic analysis of this type of problem 1s 
common and the expected dynamics are supported by observational data. With y=0, 


the jet 1s westward and the eddy rotates as before, anticyclonically. Figure 2.13 shows 
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Figure 2.9 Experiment # 3 (Jet Moved Farther Offshore) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at oy 15 yom}, (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m 
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Figure 2.10 Experiment # 4 (Topography Removed) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at dav 15 am} (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.11 Experiment # 5 (Increased Lower Layer Velocity) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at dav 15 yam, (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.12 Experiment # 6 (Lower Layer Flow Direction Reversed) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at dav 15 
interface at day 15 (m), and (d) interface at day 3 


om), (b) surface at day 30 (cm), (c) 
0 (m 
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the output from this simulation. It is readily apparent that there are some type of 
instabilities in the flow but they are substantially reduced from the cases where there 
was meridional flow. This makes sense from a linear instability standpoint and will be 
clarified in the next chapter. 

h. Experiment No. 8 (Zonal Flow Eastward) 

As with the previous experiment, the basin was not rotated into its 
standard CCS orientation. To examine the model response for a jet flowing zonally 
eastward, the basin was rotated 180°. As mentioned before, the dynamics of this type 
of configuration are well known and should help verify the model responses. Figure 
2.14 shows the output fields for this run. Of note is the fact that the flow is almost 
completely stable throughout the model run. Again, this makes sense from a 
dynamical standpoint and will also be addressed in the following chapter. 

1. Experiment No. 9 (Eddy and Topography Removed) 

In experiment # 2 the eddy was removed and in experiment # 4 
topography was eliminated. This does not totally resolve their involvement with the 
wave structure development. A requirement for baroclinic instability is a perturbation 
of the mean flow. It 1s conceivable that the eddy and the topography could be causing 
the jet to be deflected so as to undergo perturbation growth by instability theory. To 
see whether this is true, both the eddy and topography were removed for this run. As 
can be seen in Figure 2.15, the removal of both the eddy and topography did not 
hinder at all the development of the wave structures. There are some slight differences 
from the reference state run, but they are so small that it would be difficult to 
speculate as to their exact connection with the altered conditions of this run vice the 
reference one. 

j. Experiment No. 10 (Northwest to Southeast Flow) 

As can be seen in Figure 2.1, the coastline south of Cape Mendocino 
changes orientation from north-south to northwest-southeast. To examine the effects 
of this orientation on the flow the basin was rotated 120° as shown in the figure by the 
dotted box. This new configuration is important to the degree of applicability the 
model results will have to the CCS. It is obvious that the coastal jet cannot continue 
straight south after Cape Mendocino and if wave structures are not observed with the 
modeled jet going southeast it is doubtful whether the meridional flow results are true 
representations of the dynamics of the CCS. Figure 2.16 shows the model response to 


the basin rotation of 120°. The eddy was removed from this run so the wave structure 
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Figure 2.13 Experiment # 7 (Zonal Flow Westward) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at day 15 (cm), (b) surface at day 30 (cm), (c 
fnteriace at day 15 (m), or fa) interface at oy 56 (0) a ee 
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(c) (d) 
Figure 2.14 Experiment # 8 (Zonal Flow Eastward) 


Model Output Height Anomalies. 


Height anomaly fields for (a) surface at day 15 any (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.15 Experiment # 9 (Eddy and Topography Removed) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at dav 15 om (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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changes could be better seen. There are substantial differences from the reference but 
they are clearly in terms of the wave structure direction and wavelength rather than 
representing a totally different type of phenomenon. In the next chapter, the reasons 
for these structural differences will be examined. 

k. Experiment No. 11 (V, Reduced to 0.2 m]s) 

Up until now all of the model runs have used an upper laver velocity of 0.4 
m,s. This may be representative of the CCS at certain times but by no means is the 
CCS consistent. As with the previous experiment, applicability to the real world 
demands that the model parameters be as realistic as possible in both value and range. 
The velocity is reduced to 0.2 m/s thus effectively halving the vertical shear. Figure 
2.17 shows the results of this alteration. The wave structures are definitely present in a 
form similar to the reference although their strength is considerably less. This makes 
sense from an instability standpoint as the shear and thus available potential energy of 
the system has been reduced. This will be discussed in the next chapter. 

1, Experiment No. 12 (Vy; Reduced to 0.1 ms) 

As a follow on to the previous experiment, the upper layer velocity is 
reduced again by fifty percent. The results are indicated in Figure 2.18. There is some 
very slight wave development but it is considerably less than either the previous case or 
the reference state. This would intuitively mean that the reduction of the upper laver 
velocity 1s causing the system to approach a point of stability. The fact that the 
system is becoming stable with a shear still present will be covered also in the next 


chapter. 
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Figure 2.16 Experiment #10 (Northwest to Southeast Flow) 


Model Output Height Anomalies. 
Mm Cojesuiiace atedayesd cm), (Cc) 


(my 


Height anomaly fields for (a) surface at dav 15 (cm 
interface at day 15 (m), and (d) interface at day 30 
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Figure 2.17 Experiment #11 (V,; Reduced to 0.2 m/s) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at day 15 yom, (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Figure 2.18 Experiment #12 (V,, Reduced to 0.1 m/s) 
Model Output Height Anomalies. 


Height anomaly fields for (a) surface at day 15 yom) (b) surface at day 30 (cm), (c) 
interface at day 15 (m), and (d) interface at day 30 (m). 
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Il. BAROCLINIC ANALYSIS 


A. SETTING UP THE LINEAR INSTABILITY PROBLEM 
1. Instability of a Non-Zonal Flow 

Most examinations of baroclinic instability start with the assumption of a 
Zonal flow. This assumption cannot be made when dealing with the CCS where the 
basic state is definitely non-zonal and in some regions completely meridional. To 
develop an analytic model for the CCS, the quasigeostrophic equations of motron will 
be emploved. This represents a divergence from the model dynamics which are 
operating under the primitive cquations. To evaluate whether the quasigeostrophic 
equations are reasonable approximations, the Rossby number, Ro, can be calculated. 
Using U=0.4 m/s, f=0.94 x 107 st and L=Ry= 18.4 km yields Ro= 0.23 iis 
indicates that nonlinearities are appreciable but the quasigeostrophic equations should 
be within the required accuracy. 

A non-zonal flow automatically 1mphes the existence of an external forcing 
field (Pedloskv, 1979). In this case, it ts assumed that the external forcing involved 
with the CCS is wind stress curl. For a two laver system without bottom friction the 
governing quasigeostrophic potential vorticity equations in non-dimensional form from 
Pedlosky (1979) are: 


O wd CW 





See eo, en tae! 
ae 7 Ox dy Gy Ox 2 
= TG | k © curiae Be) 
Pp yf UD ec 
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‘ae * Ox dy 7 dy ax |” Y27*2 a (3.2) 


Notation is defined in the Appendix. 
Following the derivation by Pedlosky (1979) for the phase speed of a 


baroclinic wave in a non-zonal flow, consider the following basic state: 


i T = ee (3.3) 


d() 


vw, = WV, =0-- (3.4) 


The parameters u, and v, are constant velocities in the upper layer in the x and y 
direction respectively. Note that the baste state is characterized by a uniform flow 
which does not vary tn the cross stream direction. This will inherently climinate 
barotropic instability from the analytic development but duc to the length scales of the 
instabilities observed in the model runs (L> > Rq). barotropic tnstabiltty is not likely 
to be a dominant force (Kamenkovich er al., 1986). There are no boundary constraints 
on the basic flow elinunating shear induced vorticity which is present in the numerical 
model. Again, it ts felt that thrs will primarily affect barotropic instability modes. The 
wall effects on baroclinic instability are considered negligible. As Mysak et al. (1981) 


poimted out, 


. inclusion of vertical side walls at the edges of a coastal current flowing along 
a continental slope docs not significantly affect the growth rate, frequency or 
wavelength of the maximum unstable wave.... 


The real motivation for making these assumptions is that the analytic modeling of a 
non-uniform, non-zonal flow with boundary effects would be far beyond the scope of 
this presentation and would probably not yield significantly different results. Nfodel 
runs with small flows in the lower layer in both a north and south direction were so 
sinular that the assuinption of no flow tm the lower layer should be consistent with the 
numerical model. 

To examine the stability of the basic state, consider perturbations of the 


following form: 


ee ee OGL) (2y5) 
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pee (O(a ,t) (226) 


Substitution into 3.1 and 3.2 and linearization ytelds 
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Note that the wind stress curl forcing does not appear in the perturbation cquations. 
NS Pedlosky4{ 1979) describes 1m 


.. the stability of basic flows which are not solutions of the wforced equations 
of motion may be consistently considered in the context of the unforced 
perturbation cquations, without the need to consider explicitly the forces required 
to produce a basic state. 
This provides the basis for not specifying the form of the wind Stress curlameome 
Using the normal mode approach, plane wave solutions of the following form 
are assumed: 


We = A eXIGx + S260)” Simi lez (3.9) 


Substitution into 3.7 and 3.8 leads to 


AL (o-u,k-v,€)(h2+8,) + BR 


+ E(u ktv,€)] = A,F,(o-u,k-v,C) ae (3.10) 


A,LO(K*+E,) + Pk - F,(u,k+v,l)] - AF,o = 0 (3.11) 


where K* = k* + (? 


A coordinate system is now defined as shown in Figure 3.1. The angle the 
perturbation wave vector makes with the x-axis is 0 while @ represents the angle the 


shear flow makes with the x-axts. 
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Figure 3.1 Coordinate System Orientation. 


Based on the flow directions of the CCS within the numerical model domain noted in 
section II].A, @ can take on only one of two values, 270° or 300°. The coordinate 
system allows the following transformations and the introduction of the shear velocity, 


ie: 


Bk/K = B cos @ (312) 
(uk+v,0)/K = V.cos(a-0) (3.13) 
E 2 Dalz 
Where Ve (u,-+v,")"' 
The perturbation wave phase speed c 1s allowed to take on complex valucs, 
= G/K =e Sic: (3.14) 


After dividing 3.10 and 3.11 by K, use of 3.12, 3.13 and 3.14 results in, 
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AL (c-Vecos(a-0))(K7+F,) + Peosd + F,V.cos(«-8)] (alee 


- AE, (c-Vecos(u-0)) = 0 


Plea Ges) + Pcos® - F,V.cos(a-0)] - A F,c = 0 (3.16) 


Nontrivial solutions for A, and A, exist only if the determinant of their coefficients in 
3.15 and 3.16 1s equal to zero. Since the basic intent of this analytic development is 
specifically to support analysis of the model output and the CCS it is more practical at 
this point to switch to dimensional form. Using the notation of Kamenkovich er al. 
(1986), the zero determinant condition leads to the phase speed of the perturbation 


wave (Equation 3.17). 


VK*( K*+2A,)cos( a-8) - B( 2K°+A,+A,)cos0 





oe — 
2K2( K*+A,+A,) 
( (A, +A,)*B*cos*® + 2BVK*( A,-A,)cos@cos(a-8) (Gun 
- Kiv.7( 4A,A,-K*)cos?(a-8)) 17 
in 2 





2h°(K*+A,+A,) 


Where A, =F,/L*, A,=F,/L?, B=B, and all velocities from before multiplied by 
scaling velocity U. 
2. Critical Velocity and Growth Rate 
When the radicand in 3.17 is equal to zero, the wave is at the transition point 
between stability and instability. At this point, V, can be solved for. This results in an 


expression for the critical shear velocity for instability, V... 


BK7(A,-A,)cos® + 2B(A,A,)1/2[ (A,+A,)?-K*] ?74cos0 


tee K2( 4A ,A,-K*)cos(a-8) (3.18) 


-i4 


Utilizing the two basic state flow directions to match the CCS (a=270°, 
300°) and using model values of A, =2.94% 10°? m™*, A,=1.02x 10719 m‘? and 
B=2.0 x 10°F! mls! leaves V., a function of only K and 9. Figures 3.2 and 3.3 show 
the marginal stability curves for a of 270° and 300° respectively for various 
perturbation wave directions, 8. Note that curves are not shown for cases where the 
wave direction 1s perpendicular to the flow direction since by 3.18 this results in infinite 
shear velocities and no instability is possible. The form of the marginal stability is 
determined by the quadrant that 9 is in. The unstable region is shaded and represents 
the area where the shear velocity is greater than required for instability. The excess 
shear then becomes the energy source to provide a positive growth rate in those 
regions. The full implications of a meridional flow and a meridional wave direction are 
indicated by Figures 3.2b and 3.2e.. The required shear flow goes to zero and thus any 
non-zero shear will result in instability as long as the wavenumber is to the left of the 


vertical asymptote corresponding to the condition, 
K < (4A,A,)1/4 (3.19) 


which comes directly from Equation 3.18. 

With the shear conditions for instability defined, the next step is to find the 
magnitude of the growth rate in the unstable regions. Growth rate values will provide 
the necessary information to determine the characteristics of the fastest growing wave. 
For an unstable wave to grow, the radicand of 3.17 must be negative and the growth 


rate, Ke: becomes, 


[ K*v7( 4A,A,-K*)cos?( a-8) ' 
4 - 2BVK"(A,-A,)cosOcos(a-8) - (A,+A,)?82cos?9) 1/2 
a = _ 


1 





2K( K*+A,+A,) (3.20) 


As before, model values for Ay A, and B along with flow angle a reduces the growth 
rate dependence to just V,, @ and K. Figures 3.4 and 3.5 show the growth rates for @ 
= 270° and 300° respectively for velocity shears of 0.4, 0.3, 0.2 and 0.1 m/s. As 
expected, the growth rate magnitude is a strong function of the shear velocity since 


that provides the energy input into the growing wave. The maximum growth rate, 
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Figure 3.2 Marginal Stability Curves for a= 270°. 


Curves for 8 of (a) 45°, (b) 90°, (c) 135°, (d) 225°, (e) 270° and (f) 315°. 
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Figure 3.3 Marginal Stability Curves for a = 300° 


Curves for 8 of (a) 75°, (b) 120° 


(c) 165°, (d) 255°, (e) 300° and (f) 345° 
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Figure 3.4 Growth Rates fora 


Growth rates ( Xx 107° sh) for velocity shears of (a) 0.4, (b) 0.3, (c) 0.2 and (d) 0.1 m/s. 
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which would correspond to the fastest growing. wave, occurs at different K-0 
intersections in each different shear flow case. This means the shear flow not only 
affects the magnitude of the growth rate but also the direction and wavenumber of the 
perturbation wave. As a result, any variation in the external forcing which is 
producing the shear flow would cause major alterations to the fastest growing wave. 
Another interesting point is that the growth rates for @=270° are higher in all shear 
conditions than for a=300° although the 6-K plane positive growth rate areal extent 
is smaller. From the marginal stability curves, it is apparent that a meridional flow 
could have perturbations for which any non-zero shear would support instability. As a 
result, the additional claim that meridional flow yields the largest growth rate for a 


given shear can be made. 


B. APPLICATION OF ANALYTIC RESULTS TO MODEL OUTPUT 

The instability problem has now been formulated to approximate the numerically 
modeled situation. The next step is to examine the instabilities present in the model 
output fields as noted in section II.B.2 in terms of the previous section’s results. If 
there are no inconsistencies, then the instabilities are most likely due to baroclinic 
instabilitv of the basic flow. 

From section II.B.2, it is obvious that there is a similarity of the output fields 
With respect to apparent wave-like structures. Not only do the same features appear in 
most of the runs, but their shape, direction and magnitudes are also similar. Table 3 
shows the average values of relevant parameters from the model runs (where observed) 
for days 10, 20 and 30 for flow direction of a=270°. Table 4 shows the same 
parameters for flow direction of a= 300°. The wave data for both tables was taken 


from model runs where the shear flow was close to 0.4 m/s (+ 0.01 m/s). 


TABLE 3 
WAVE PARAMETERS FROM MODEL OUTPUT FIELDS (a= 270”) 


DAY Q (deg.) K (m’}) L (km) 


10 os 2.0 x 10° 310.0 
249 2.6 x 10° 240.0 
255 3.1 x 10° 200.0 
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TABLE 4 
WAVE PARAMETERS FROM MODEL OUTPUT FIELDS (a= 300°) 


DAY Q (deg.) K (m‘!) L (km) 


10 238 3 105° 275.0 
20 260 3.0 x 107° 209.0 
30 269 3.2 x 107° 198.0 





These days were selected because they represent the appearance (day 10), the 
growth (day 20) and the maturity (day 30) of the waves. By mapping the wave data 
from Tables 3 and 4 onto instability plots from section III.A.2, an evaluation of their 
relationship to the analytically derived baroclinic wave characteristics can be made. 

Marginal stability curves are shown in Figure 3.6 for the flow angles and wave 
angles noted in Tables 3 and 4. The wavenumbers for those angle values are indicated 
and marked by the day the data is associated with at its intersection with the V,.=0.4 
ms line. Of note is the fact that all of the plots show the waves within the unstable 
region for baroclinic instability. Additionally, as the waves shift in wavenumber 
towards shorter wavelengths, their shear excess increases to its maximum where the 
wave meets the asymptote corresponding to Equation 3.19. 

Growth rate plots are shown in Figure 3.7 for velocity shears of 0.4 m/s. The 
parameter positions corresponding to davs 10, 20 and 30 from Tables 3 and 4 are 
indicated to show the path the waves took as they developed. From Figure 3.7 it is 
apparent that the model output waves fall in the regions of positive growth rate but do 
not fall on or near the peak growth rate which would correspond theoretically to the 
position of the fastest growing wave. Also, the growth rate 1s low at day 10, reaches a 
maximum at day 20 and 1s decreasing by day 30. It must be kept in mind that these 
growth rates represent the analytic values given the wavelength and wave direction of 
the model output waves. They are not necessarily the same as the growth rates of the 
actual waves in the model output. The basic effect of the growth rate is to increase the 


amplitude of the growing wave in the following manner: 


A = Agel Keid (3.21) 
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Figure 3.6 Marginal Stability Curves for Model Output Wave Characteristics. 
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Figure 3.7 Growth Rate Positions of Model Output Waves for Shear Velocity of 0.4 m/s. 


Parameter data from (a) Table 3 and (b) Table 4 ploted on growth rate ( x 10°9 sly 
fields from Figures 3.4.a and 3.5.a. 
where A is the time varying amplitude and Ap is the initial amplitude. Rearranging 


Equation 3.21, 
En(A - Ap) = ket (2) 


yields an expression of the form, y= mx where m represents the slope (kc;). By plotting 
the natural log of the amplitude difference from the initial versus time, the growth rate 
can then be determined by the slope of the curve. Figure 3.8 shows the amplitude 
changes of the reference state waves during the model run. 

Immediately apparent is the fact that the growth rate varies with time. It 1s low 


at first and grows to a maximum and then decreases. The information from Figure 3.8 
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Figure 3.8 Reference State Wave Amplitude Growth. 


can be easily converted into a plot of growth rate versus time. This is shown in Figure 
3.9. The temporal variation of the growth rate is more clear and shows similar form to 
the growth rate progression observed in Figure 3.7.a. 

Experiments 7 and 8 involved zonal flows. These were not meant to simulate the 
CCS but to evaluate the model response against theory. For a zonal flow, the wave 
growth must offset one hundred percent of the stabilizing planetary vorticity gradient, 
whereas for a meridional flow the gradient’s effects are zero due to the Pcos@ term. 
Figure 3.10 shows the marginal stability curves for an eastward and westward zonal jet. 
The curves are not for any specific 8 because the 8 dependence cancels out of equation 
3.18 for a zonal flow. Immediately apparent is the fact that the easterly jet should be 
more stable than the westerly jet since it has a much higher munimum shear velocity. 
This agrees with the experimental results where the easterly jet remained stable through 


the entire run and the westerly jet developed slight instabilities. Based on the curves 
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Figure 3.9 Growth Rate Variations of Reference State Waves. 


for the zonal flow coupled with the fact that the shear velocity was 0.4 m/s, one might 
expect instability wave growth for both zonal flows similar to the reference run 
although this was not observed. The increase in the required V,, for the zonal flows 
due to the B effect cuts down on the excess shear available for growth and apparently 
this excess is not quite great enough for wave growth. This points to frictional effects 


and will be discussed in the next section. 


C. FURTHER ASPECTS OF THE INSTABILITY PROBLEM 

The assumptions made to support the application of baroclinic instability theory 
to the numerical model results were necessary for any quantitative analysis to be made. 
Generally, there was very good agreement between linear instability theory and the 
model output which supports the contention that baroclinic instability is in fact the 
dominant dynamical process in the model and in the CCS. Nevertheless, there were 
discrepancies between the model siinulations and the analytic predictions and these 


need to be addressed. 
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Figure 3.10 Marginal Stability Curves for Zonal Flow. 


Curves are for @=270° and flows to (a) east and (b) west. 


1. Effects of Friction 

Without friction, the vertical phase shift of a baroclinic wave opposite the 
vertical velocity shear 1s an indication that Cre Q and the wave should grow (Pedlosky, 
1979). However, when frictional forces are introduced, the wave can take on the same 
phase shift opposite the velocity shear and still be neutral. Instead of supplying energy 
for wave growth, the extracted energy is spent offsetting dissipation. Thus, in order to 
grow in the presence of friction, the baroclinic wave must have a greater phase sluft 
than would be the case without friction. Since there is a maximum phase shift 
associated with maximum energy extraction, the growth rate of a baroclinic wave with 
friction should be less than that for a wave without friction. 

In Figures 3.4 and 3.5 it was observed that changes in the shear velocity 
resulted in maximum growth rate value and 9-K positional changes. Frictional forces 


would tend to reduce the available energy for growth which is equivalent to a reduction 
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in the shear velocity. For this reason, the numerical output waves could not be 
expected to fall on the maximum growth rate parametric positions predicted by the 
non-frictional analytic development. Another complicating factor is the way in which 
friction was parameterized in the numerical model. By using a Ay,V?V form, friction 1s 
a strong function of the velocity field. Since the velocity field in the domain changed 
drastically as the baroclinic waves developed, the frictional forcing must also have 
changed drastically. Thus, the frictional alteration of the maximum growth rate 0-K 
position would take on temporal variability. This might explain one reason the model 
output waves were observed to vary in wavenumber and direction through the model 
runs. 
2. Non-linear Instability Effects 

In linear instability theory the perturbation wave is unbounded as long as it 
meets the necessary conditions for instability which are defined in terms of the basic 
state. Obviously, this is not realistic as the waves grow to finite amplitude. Non-linear 
theory explains how advection affects the wave growth. 

As Pedlosky (1979) pointed out, it is the basic state in the absence of the 
perturbation wave that must be used as criterion for instability vice any instantaneous 
flow profile. Variations in the flow profile as the wave develops, however, does have 
an effect on the wave. When the wave reaches a point where the instantaneous mean 
flow would not support instability, the wave growth rate is at its maximum and will 
Start to decline. Although declining, the growth rate 1s still positive and the wave is 
still growing even though the instantaneous mean flow might support no 
instability. The change in growth rate is accompanied by a reduction of the phase shift 
of the wave in the vertical until growth is stopped as the phase shift reaches zero. The 
process does not stop here but continues into the phase orientation that causes wave 
decline. As with the growth, the decline wili hit a maximum and reverse and the whole 
process of growth and decay will oscillate back and forth (Pedlosky, 1979). Thus, the 
variations in the wave growth rate observed in Figure 3.9 and inferred from Figure 3.7 
could well be indicative of non-linear effects present in the model and presumably, the 
WES. 

When considering the effects of friction, variations in the available energy to 
the wave were correlated to variations in the parametric characteristics of the 
maximum growth wave. Non-linear theory predicts similar variations in the available 


energy due to the changes in the mean flow brought about by the perturbation wave. 


> 


Hence, it is another possible causative factor for the. model output, since the model is 
non-linear, to exhibit variability in baroclinic wave wavenumber and direction over 
time. 

3. Mixed Instability Effects 

Most treatments of the instability problem are either focused on baroclinic or 
barotropic instabilitv. In the CCS and the numerical model, flows have both 
horizontal and vertical shears. [t was assumed in the analytic development of the 
baroclinic instability problem in section III.A.1 that barotropic instability was not the 
dominant instability process based on the scales of the observed disturbances. 
However, it may be that barotropic instability does have a slight influence on the 
baroclinically unstable waves. 

Pedlosky (1979) showed that for a system with vertical and horizontal velocity 
shears, the characteristics of the most unstable wave can only be derived with 
exhaustive calculation given a detailed spatial distribution of velocity and potential 
vorticity. Changing the parameters of the flow can have a strong inpact on the type of 
instability present and the wavelength. It is believed that throughout the numerical 
simulations, baroclinic instability remained the dominant instability process but 
barotropic instability may have altered the characteristics of the fastest growing wave. 

4. Baroclinic Instability with Complex Wavenumbers 

In developing the baroclinic instability problem, a particular perturbation 
form was assumed. This form is common to nearly all studies of this type of instability 
and is shown by Equations 3.9 and 3.14. The crucial point is the selection of a real 
wavenumber and a complex phase speed. 

Hogg (1976) examined analytically and numerically the baroclinic instability 
process in a zonal flow where the perturbation form allowed both complex phase speed 
and complex wavenumber. In general, his results indicated a wavenumber cutoff 
between spatially growing (k;>0) and temporally growing (c,>0) waves. Temporally 
growing waves were confined to long wavelengths while spatially growing waves were 
confined to the shorter wavelengths. 

In a study of Gulf Stream meanders, Watts and Johns (1982) used statistical 
analvsis of cross stream data and derived values for both temporal and spatial growth 
rates of the meanders. Contrary to Hogg (1976), Watts and Johns’ (1982) results 
indicated that spatial and temporal growth rates can occur in similar wavelength 
regimes. Figure 3.11 shows the growth rate versus wavenumber relation for both 


spatial and temporal growth found by Watts and Johns (1982). 
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Figure 3.11 Spatial and Temporal Wave Growth of Gulf Stream Meanders 
(Watts and Johns, 1982). 


There is evidence both analytically and observationally that spatial growth 
may be involved in the baroclinically unstable wave. Although Ilogg (1976) 
concentrated on spatial growth (k;>0) and Watts and Johns (1982) observed only 
temporal and spatial growth (c;>0, k;>0), there is nothing in the theory that prevents 
a combination of temporal growth and spatial decline (c;>0, k;<0). This might also 
explain the variation in the observed wavelengths in the model output from long to 


short while the waves were growing in time. 
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IV. DISCUSSION AND CONCLUSIONS 


A. ANALYSIS OF RESULTS 
1. Comparison of Results with Observational Data 

As mentioned in the introduction, filament structure, spatial distributions and 
temporal variabilities are well Known. [n order to validate the results of this study, 
there should be general agreement with the existing information and if discrepancies 
are present there should be some theoretical basis for the differences. In Table 5 
values are shown for the range and magnitude of the observed versus the modeled 
characteristics. The value for the observed wavelength and direction is from Ikeda and 
Emery (1984) and the observed offshore extent 1s from Brink and Hartwig (1985). The 


rest of the observed values are from Flament er al. (1985). 










TABLE 5 
| FILAMENT CHARACTERISTICS, OBSERVED VS WiIGBErER 


PARAMETER OBSERVED MODELED 
Wavelength 100-400 km 198-310 km 
Wave Direction ~ 270° 215-269° 

Off-shore Extent 100-500 km 100-800 km 







Maximum Off-shore Velocity =~'0) Ss O25 mys 
Maximum On-shore Velocity ~0.35 m/s 0.4 m/s 


Growth rate (e-folding time) ~ 1-2 weeks — 10nd 






It would seem that the only parameter that deviated significantly was the wave 
direction. The reason for this will be discussed in the next section. Otherwise the 
agreement strongly supports the contention that the modeled phenomena are in fact 
representative of filaments in the CCS. 
2. The Filament Formation Process 
The filament formation process is quite complex. Although baroclinic 
instabilitv is clearly the dominant mechanism, there are other forces at work that 


cannot be neglected. 
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In the first phases of growth, the wave is-affected by the east-west horizontal 
velocity shear. Close to the jet, the wave front is advected strongly to the south. 
Farther offshore, the jet effect on the wave propagation decreases due to the mean 
velocity gradient. This causes the wave vector direction to be deflected from the 
seemingly more favorable southerly direction to a near westerly direction. As the wave 
gains energy from the vertical shear of the mean flow. advection of the mean flow ina 
zonal direction starts to break down the east-west horizontal velocity shear and the 
wave direction rotates towards its more optimum direction and its wavelength shortens 
to maximize the growth rate. Non-linearities become dominant as the wave reaches 
such a magnitude that the instantaneous mean flow 1s stable and the wave growth 
starts to decline. This 1s as far as the waves were modeled, but it would be expected 
that the waves would go into an oscillatory growth-decline mode until damped by 
friction. 

Obviously, this whole analysis 1s contingent on a constant wind stress forcing 
for periods on the order of months which 1s hardly realistic. It can be expected that 
constant forcing mav persist for weeks at a time in which case the initial growth and 
development holds true. 

From an observational standpoint, the filaments observed are nearly always 
oriented perpendicular to the flow indicating a southerly wave direction. This does not 
mean that the described theory of development is necessarilv flawed, it just means that 
by the time the waves are observable, they have probably already twisted around into 


their more optimum configuration. 


B. IMPORTANCE OF RESULTS 

There have been numerous studies of zonal instability processes and a few of 
meridional situations but there is a large and distinct scarcity of studies of flows in- 
between. When one considers currents such as the Gulf Stream, the Kuroshio and, of 
course, the CCS which have flows that are not always east-west or north-south, it 
becomes evident that more attention should be placed on this discrepancy. 

As demonstrated, the behavior of a flow that is non-zonal is dramatically 
different from a zonal flow. In a zonal flow, the energy releasing horizontal 
temperature gradient and the stabilizing planetary vorticity gradient are always working 
against each other. In a meridional flow, energy release can take place without any 


component of § acting on it; hence, any shear above the dissipation level may produce 
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instability. For the in-between case, the wave must assume a balance between the 
maximization of energy release and the dampening effects of the planetary vorticity 


gradient. 


C. RECOMMENDATIONS FOR FUTURE STUDIES 

There are several aspects of this modeling work that could be modified to better 
simulate the CCS. A better resolution model would help significantly to bring the 
barotropic instability process into view. This is particularly important to the filament 
situation since observational evidence has shown instabilities along the filament frontal 
region of the scale of Rg that are claimed to be the result of barotropic instability 
(Flament er al., 1985). A model with more layers would give the vertical resolution 
necessarv to properly simulate the CCS in all of its configurations, particularly its late 
summer profile of a strong southward surface jet and a poleward undercurrent all 
above a relatively quiescent deep water region. The addition of some temporal 
variability to the coastal jet strength would help to better understand wind stress 
variation effects on existing filaments and the growth of future filaments. Continued 
intensive in situ and satellite gathering of oceanographic data is vitally important to 
understanding the physical processes necessary to better reproduce the observed 
filament structures numerically. To thoroughly determine if baroclinic instability is 
occurring in the CCS, an elaborate experiment would have to be designed to resolve 
spatially and temporally the mesoscale motions. An experiment comparable to the 


POLY MODE experiment in the North Atlantic would be necessary. 


D. CONCLUSIONS 

Baroclinic instability is a prime generating force for filaments formed south of 
Cape Mendocino in the CCS. Barotropic instability along with non-linearities, 
frictional effects and possible spatial wave growth may help to modify the 
characteristics of the growing baroclinically unstable waves. 

Filament growth is strongly dependent on the vertical velocity shear, the vertical 
structure and the mean flow direction. Any attempt to model filaments must carefully 
parameterize these quantities if accurate and thus useable results are to be obtained. 
Indeed, for any numerical simulation of currents with unknown instabilities, their 
values are crucial to correct simulations. Improper understanding of their spatial and 


temporal characteristics can seriously affect the validity of the model results. 
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APPENDIX 
SYMBOLS AND NOTATION 


Laplacian lateral friction coefficient 

Amplitude coefficient for normal mode solution for nth layer 
Phase speed of perturbation wave 

Variable depth of topography 

Upper layer thickness 


Froude number for i'” layer nae eel: 
Coriolis parameter for mean latitude 
Gravitational acceleration 

Reduced gravitational acceleration = 2(P,-P)/P; 
Instantaneous laver thickness 

Upper (:= 1) and lower (1:= 2) layer mean thickness 

Unit vector in the z direction 

Perturbation wave number in the x direction 

Perturbation wave number in the y direction 


k? +0? 


Total wave number for perturbation wave 
Horizontal scale length 

Pressure in the upper layer = g(h, +h, +d) 

=p,-gh, 

=(¢'H,H,/f,°(H, +H,)]!/? 
=U/fL 


Pressure in the lower layer 

First internal Rossby radius of deformation 
Rossby number 

Scale velocity 

Velocities in the x and y directions 
Transport in the x and y directions 

Shear Velocity = (un ate a ei 
Critical shear velocity for instability 

Cartesian coordinates directed N and W respectively 

Angle mean flow vector makes with x axis 

= df/dy 


=§,L7/U 


Variation of Coriolis parameter with latitude 
Scaled variation of Coriolis parameter with latitude 
Basin rotation angle from east 


Grid spatial resolution 
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Model time increment 


Angle perturbation wave vector makes with x axis 
Dimensionalized Froude number for i layer 
Density in jth layer 

Mean density 

Frequency of perturbation wave 

Wind stress vector 

Perturbation stream function for it layer 

Basic state stream function for ith laver 

Variable stream function for it? laver 

Gradient operator 


Laplacian operator 
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= P12 
= FIL 


= 0/0x + d/dy 
= G*/6x? + 6*/dy? 
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